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Abstract 

On the basis of linear programming, new sets of entanglement witnesses (EWs) for 3 <8> 3 
and 4 ® 4 systems are constructed. In both cases, the constructed EWs correspond to the hyper- 
planes contacting, without intersecting, the related feasible regions at line segments and restricted 
planes respectively Due to the special property of the contacting area between the hyper-planes 
and the feasible regions, the corresponding hyper-planes can be turned around the contacting area 
throughout a bounded interval and hence create an infinite number of EWs. As these EWs are 
able to detect entanglement of some PPT states, they are non-decomposable (nd-EWs). 

Keywords: Entanglement Witness, Feasible Region, PPT states. 
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Entanglement which appears only in composite quantum systems, is the main difference between 
quantum mechanics and classical physics. It is used as a physical resource to realize various quantum 
information and quantum computation tasks such as quantum cryptography, teleportation, dense 
coding, and key distribution 03 El El. The fundamental problem of entanglement theory is the finding 
out separating boundary of the entangled states and separable ones. The celebrated Peres-Horodecki 
criterion based on positive partial transpose (PPT) determines this boundary for 2 <g) 2 and 2 <g) 3 
cases @] |5l but it has no efficiency for PPT entangled states appearing in the higher dimensional 
systems. On the other hand, distinguishing between the PPT entangled states which can not be 
distilled and those which can be distilled is an important problem in quantum communication field 
lEfl . The most general approach to study the entanglement of quantum states in higher dimensional 
physical systems is based on the notion of entanglement witnesses (EWs) Q [SJ. The EWs are 
essential tools in entanglement theory since it has been shown that for any entangled state there 
exists at least one EW which detects its entanglement EUl. A Hermitian operator W is said to be 
an EW if and only if for all separable states p sep , Tr(W p sep ) > and at least for one entangled state 
p ent , Tr(W p en t) < ( one says that p en t is detected by W). Clearly, the construction of EWs is a 
hard task. Although it is easy to construct a Hermitian operator W which has negative expectation 
value with some entangled states but it is very difficult to check that its expectation values with all 
separable states are non-negative. So, several approaches for constructing EWs have been proposed. 
It has been shown that among these approaches, the linear programming, as a special case of convex 
optimization [U0l[TT| . is a very useful one [|T2l IT3l IT4l IT51I . In linear programming method, the main 
problem is the determination of boundaries of a convex set known as feasible region. The feasible 
region is associated to the convex set of separable states. The boundaries of the feasible region 
are characterized by the operators used in constructing the EWs. When the feasible region was 
specified then every hyperplane tangent to the feasible region corresponds to an EW. A hyperplane 
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may be contact the feasible region at a point, line segment, restricted plane, etc. If the boundary is 
differentiable on the contact area, the hyper-plane has to be fixed, otherwise it can be turned around 
the contact area throughout a bounded interval. The EWs constructed so far via linear programming, 
correspond to the fixed tangent hyper-planes. When a boundary of the fesible region is a hyper-plane 
by itself, it is called an exact boundary. 

In this paper, by choosing suitable operators, new sets of EWs are constructed for 3 <g> 3 analyti- 
cally and for 4 <g> 4 numerically. These EWs correspond to the hyper-planes which can turned around 
the contact line segments or restricted planes for the 3 <8> 3 and 4 <g) 4 cases respectively. It is shown 
that these hyper-planes can be turned around in a bounded interval such that they do not intersect the 
feasible region. In this way, we obtain an infinite number of hyper-planes which are in contacting 
with the feasible region and hence an infinite number of EWs for 3 <g> 3 and 4 <g> 4 systems. Finally, we 
show that these EWs are nd-EWs since they are able to detect the entanglement of some PPT sates. 

The paper is organized as follow: In section 2, for 3 ® 3 case, we obtain a part of the boundary 
of the feasible region as a exact boundary. The intersection of this boundary by the other boundaries 
is three line segments. We see that the planes which contact the feasible region only at these line 
segments, correspond to new EWs. Section 3 is devoted to extend the approach of section 2, to the 
4 <g) 4 systems. It is ended by a brief conclusion. 

2 Entanglement Witnesses For 3 ® 3 Systems 

In this section we are going to introduce a method based on linear programming to construct new 
EWs on 3 <S> 3 Hilbert space which are able to detect the entanglement of some PPT states. When 
one deals with linear programming, the aim is optimization of a linear function under some linear 
constraints. The linear constraints define a region called feasible region. So to reduce the construction 
of EWs to a linear programming problem, the related feasible region must be characterized. As 
mentioned in introduction, the shape of the boundaries of the feasible region depends on the operators 
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chosen for constructing of EWs. So in this paper, we choose the following operators 

1 = |(|01)(01| + |12)(12| + |20)(20|) 

2 = |(|02)(02| + |10)(10| + |21)(21|) (2.1) 

, where 

|^) = i=(|00) + |ll) + |22)). 

It is seen that these operators are orthogonal to each other, i.e. OiOj = 0, i ^ j = 1,2,3, and 
2 = nOiIl, 3 = n0 3 n where II is an operator which permutes the particles. We denote the 
expectation values of the above operators with pure product states \a) <8> \(3) by px, p 2 , P3 respectively 
as 



Pi = !(K| 2 |/3i| 2 + |ai| 2 |/? 2 | 2 + |a 2 | 2 |A) 



2\ 



P2 = U\a \ 2 W + |«i| 2 |A>| 2 + |« 2 | 2 |A| 2 ) (2-2) 



P3 = k\a (3o + OLiPt + a 2 /3 2 \ 



2 



The pi, p 2 and p 3 can be considered as components of a point lying in the three dimensional Euclidean 
space. Hence, the set of all pure product states form a region in this Euclidean space. Since every 
separable state can be written as a convex combination of pure product states, the feasible region is 
the convex hull of the mentioned region. Now let us go to figure out a part of the boundary of the 
feasible region which has essential importance in constructing our EWs. As the maximum value of 
Px, p 2 and p 3 is |, the points (|, 0, 0), (0, |, 0) and (0, 0, |) which correspond to the product states 
listed in the table 1, are vertices of the feasible region. The plane which passes through these points 
has the following equation 

3(pi+p 2 +p 3 ) = l (2.3) 

To find out that this plane is a boundary of the feasible region, the linear function on the left hand 
side of the (12.31) must be maximized with respect to the product states. If the maximum value is one 
then the plane is a part of the boundary of the feasible region, otherwise it is not. Calculation shows 
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that the maximum value is greater than one and hence this plane is not a boundary for the feasible 
region. The product which gives the maximum value is 



i 

v/3 



(A 

1 

V 1 / 



i 

V3 



and its corresponding point is (~, ~, |). This point is another vertex of the feasible region which also 

given in the table 1. The plane tangent to the feasible region at this point is 

5 



3(pi +P2 + Ps) 



(2.4) 



This means that for all separable states, we have always 3(pi + P2 + Ps) < §■ hi the next step, we 
note that the equation of the plane passing through the vertices (|, 0, 0), (0, |, 0) and (|, |, |) is 

3(pi + P2) + Ps = 1 (2.5) 

The maximum value of the left hand side of this equation with respect to the product states is one 
so this plane is an exact boundary of the feasible region and any separable state always satisfies the 
inequality 3(pi + p-z) + p% < 1. The existence of this exact plane makes a good possibility for 
constructing a new type of EWs. The line segments surrounding the plane (|2.5I) . are also on the 
boundary of the feasible region. The sets of planes which are in contact with the feasible region at 
the mentioned line segments and not intersect the feasible region can be corresponded to new type 
of EWs which we desire to construct them. To this aim, we proceed as follows. First, we choose the 
vertices (0, ~, 0), (0, 0, |) and (|, ~, |). The plane passing through them is 

3(-pi+p 2 + p 3 ) = 1 (2.6) 
Maximizing the left hand side of this equation with respect to the product states, we obtain the 



maximum value | at the following state 
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V2 
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So the plane (12.61) is not a boundary plane and the corresponding point of the above pure product 
state, i.e. ^, j), is a vertex of the feasible region. The plane tangent to the feasible region at this 
point is 



3(-pi +P2 +Ps) = J 



(2.7) 



In the second step, we choose the vertices (0, |, 0), (|, |, |) and ^, |). The plane passing through 
them is 



3pi + 9p 2 + 5p 3 = 3 



(2.8) 



which it is not a boundary of the feasible region because its maximum value with respect to the prod- 
uct states is not 3 but ||. This maximum value is achieved at the point ^) corresponding to 



the following product state 



( o N 

V2 



\ V14 ) 

The plane tangent to the feasible region at this point is 



I o \ 

14 



3pi + 9p 2 + 5p 3 



73 
24 



(2.9) 



Finally, we choose the vertices (0, 0, |), (^, |, |) and ^, ~) and pass a plane through them. The 
plane is 



3(— pi +p 2 ) + 6p 3 = 2 



(2.10) 



Maximizing with respect to the product states shows that this plane again is not a boundary plane 



because the point ^) corresponding to the product state 



' o N 



V 



To 



( o \ 



/ 



V 
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, gives the value Hence, the tangent plane at the point (Jf, -|) is 



3(-£>i + P2) + 6p 3 



64' 192' 16' 

17 



(2.11) 
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Now we find the intersection of the planes (12.91) . (12.1 II) and the plane ps = | which is the vertex 
(i |, |). Clearly this vertex does not belong to the feasible region as it can not be obtained from 
any product state. Therefore, the equation of the plane passing through the vertices (|, ~, ~), (0, |, 0) 
and(^,|,i)is 

3pi + 6P2 + 6p 3 = 2 (2.12) 

This plane has intersection with the plane (|2.5I) throughout the line segment passing through the 
vertices (|, |, |) and (0, |, 0). In fact, the mentioned plane is tangent to the feasible region at the 
above line segment. Obviously, the plane (12.51) which is an exact boundary plane of the feasible 
region, can be rotated around the line segment and coincide to the plane (12.121) . However, the plane 
(12.121) is the final limit for the rotation of the exact plane (12.51) since if we rotate it much more, the 
resultant plane intersects the feasible region. The permutation of particles changes the plane (12.121) 
into the plane 

6p 2 + 6pi + 3p 3 = 2. (2.13) 

This plane which passes through the vertices (|, §, §), (§, 0, 0) and (|, j^, |), is tangent to the feasible 
region at the line segment passing through the vertices (|, |, ~) and (|, 0, 0). It is clear that the plane 
(12.131) can also be obtained by rotating the exact plane (12.51 ) around this line segment. By the same 
argument sketched for the plane (12.121) . the rotation is bounded by the plane (12.131) . To summarize, 
the planes obtained by rotating the exact plane around the line segment passing through the vertices 
(|, |, |) and (0, |, 0) can be written in the following parametric form 

^ + 3p 2 + (2--^)p 3 = l, (2.14) 

a 6a 

where a varies in the interval [|, |] with a = | for the exact plane (12.51) and a — | for the plane 
(12.121) . Under the permutation of particles, the parametric plane (|2.14l) is transformed to the plane 

Do 1 

^ + 3pi+ 2-— p3 = l, (2.15) 

a 6a 

which represents a set of planes obtained by rotating the exact plane (|2.5I) around the line segment 
passing through the vertices (|, |, |) and (|, 0, 0). Consequently, any separable state satisfies the 
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inequalities ^ + 6p 2 + (2 - ^)p 3 < 1 and f + 6pi + (2 - ±)p 3 < 1 for a £ [§, §]. Finally, consider 
the parametric plane 

3bi+p 2 ) + -P3 = l, (2.16) 

a 

for all a £ (— oo, 0) (J[l, °°) an d the plane p 3 = for a = 0. This parametric plane is representative 
for a set of planes which do not intersect the feasible region and obtained by rotating the exact plane 
(12.51) around the line segment passing through the vertices (|, 0, 0) and (0, |, 0). Therefore for any 
separable state, we have 3(pi + p 2 ) + ^p 3 < 1 for a ^ and p 3 > for a = 0. Now the description 
of the feasible region for our purposes is complete and hence we are going to introduce the EWs 
corresponding to the planes discussed above. To this aim, let us consider the following Hermitian 
operator 

W a = h ® h ~ -Oi - 30 2 - (2 - -L)0 3 (2.17) 
a 6a 

which is correspond to the parametric plane (12.141) . By construction, the expectation value of this 
operator with respect to the all separable states is positive hence it can be an EW for a e The 
permutation of particles transforms this EW to the following one 

W' a = I 3 ® I 3 - -0 2 - 3d - (2 - J-)O a (2.18) 
a 6a 

which is correspond to the parametric plane (12.151) . So it has positive expectation value with respect 
to the all separable states therefore it can also be an EW. Finally, consider the Hermitian operator 

W'i = h®h- 3(0! + 2 ) - -03 (2.19) 

a 

with a £ (— oo, 0) |J[1, oo) corresponding to the parametric plane (12.161) and the operator W = 3 
corresponding to the plane p 3 = 0. The operators W'^ and W are both positive so they can not serve 
as EWs. For a — |, W a = W' a is a positive operator and we claim that for a £ (|, |] the EWs W a 
and W a are nd-EWs. For this purpose, let us consider the following state 

p = a 1 1 + a 2 2 + a 3 3 (2.20) 



Entanglement Witnesses 10 

where aj > for i = 1, 2, 3 and ai + 02 + 03 = 1. It is easy to see that when a\a 2 > 03, then p is 
PPT. The expectation values of W a and W a with this state are 

Tr(W a p) = {a 1 -a 3 ){l-±-) (2.21) 

da 

and 

Tr(W' a p) = (a 2 -a 3 )(l-^) (2.22) 

So the entanglement of the PPT state (12.201) is detected by the EWs W a for < a x < a 3 < ± and 
I < a 2 < 1. It is also detected by the EWs W a for < a 2 < a 3 < | and | < eti < 1. As a special 
case of the state (|2.20l) we take the following state which was introduced in lfi~6l 

p, = yOx + + 2 -0, (2.23) 

with < p < 5. This state is separable for 2 < (3 < 3 and not PPT for /3 > 4 and < 1. So the 
equation (12.211) becomes 

Tr{W a p,) = {P-2){\-^-) (2.24) 

The W a for all a G (§, §] detects the entanglement of p p for (3 < 2. So the is PPT entangled for 
1 < f3 < 2 and free entangled for < (3 < 1. On the other hand, the equation (l2.22l) becomes 

Tr(W' aPfi ) = (3-f3)^--^-) (2.25) 

It is also seen that W a for all a G (|, |] detects the entanglement of p p for /3 > 3. So the is PPT 
entangled for 3 < (3 < 4 and free entangled for 4 < /3 < 5. Therefore W a and are nd-EWs for 
all a G (|, |] and our claim is proven. To make a comparison among the EWs W a for a G (|, |], it 
is seen from the equation (12.241) that the EW W2 , i.e. 

W2 = h ® h - \o x - 30 2 - ^0 3 , (2.26) 

has the best detection of entanglement of pp. The similar comparison can be made among the EWs 
W' a for a G (|, |]. The result is that the EW H^, i-e. 



= J 3 ® J 3 - - 30x - ^0 3 (2.27) 
3 z z 
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is the best one. This EW is the same one that obtained numerically by Doherty and et.al. in ifTTl . 
To measure the EWs W a and W' a experimentally, they should be, in fact, decomposed into a sum of 
locally measurable operators lfi~8~l [i~9~l . There exist the following decompositions for EWs W a and 

W' a as 

W a = — ^ h ®h + ~~77, — (A 2 ® A 2 - Xi ® Ai + A 5 g> A 5 - A 4 <g> A 4 + A 7 g> A 7 - A 6 <g> A 6 ) 

27a 18a 

3a — 5,. , , , . V3(3a — 1),_ , . , N 

- — — (A 3 ® A 3 + A 8 ® A 8 ) + — y — ^(A 3 <g A 8 - A 8 ® A 3 ) (2.28) 

3oa 12a 

and 

W a ' = U ^ , 2 h ®h + 6 " / (A2 ® A 2 - Ai ® Ai + A 5 ® A 5 - A 4 ® A 4 + A 7 g) A 7 - A 6 ® A 6 ) 
27a' 18a' 

3a' — 5.. , , , . v3(3a'— 1),. , . , s 

- Qfi . (A3 ® A 3 + A 8 ® A 8 K —— / — '- A 3 ® A 8 - A 8 <g> A 3 ), (2.29) 

3oa 12a 

where AjS, i — 1, 8, are the basis for the sw(3) Lie algebra [201. Therefore each EW requires ten 
measurement settings which should be measured locally by Alice and Bob simultaneously. 



3 Entanglement Witnesses For 4 ® 4 Systems 

In this section we extend the approach of the previous section for 3 <8> 3 systems to the 4 <g> 4 ones. To 
this aim, we introduce four operators as 

01 = |(|01)(01| + |12)(12| + |23)(23| + |30)(30|) 

2 = i()02> <02| + |13)(13| + |20)(20| + [31> <31|) 

(3.30) 

3 = K|03)(03| + |10)(10| + |21)(21| + |32)(32|) 
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where 



M = i(|00) + |ll) + |22) + |33)) 



Pi 


= |(l«o| 2 |/?i 


2 + 




2 |/? 2 | 2 + 


«2 


2 |/3 3 


2 + NW) 


P2 


= KKH/?2 


2 + 


\ati 


2 |/? 3 | 2 + 


«2 


2 |/?o 


2 + l« 3 | 2 |/3i| 2 ) 


P3 


= |(l«0| 2 |/?3 


2 + 


\0L\ 


2 IA)I 2 + 


«2 


2 |A 


2 + H 2 |/3 2 | 2 ) 



It is clear that 3 = nOiII, 2 = n0 2 n, A = TlOJl and 0*0, = 0, i ^ j = 1, 2, 3, 4 where n 
is the permutation operator which permute the particles. The expectation values of these operators 
with respect to the product state \a) <g) \/3) are shown by p l5 p 2 , P3 and p 4 respectively as 



(3.31) 



p A = i\a (3 + «iA + 02/^2 + a3/? 3 | 2 - 
As before, the p l5 p 2 , P3 and p 4 can be considered as the components of a point in the feasible 

region lying in the four dimensional Euclidean space. Again, we determine a part of the boundary of 

the feasible region which is suitable for our purposes. By considering the points (j, 0, 0), (0, |, 0, 0), 

(0, 0, |, 0) and (0, 0, 0, |) as the vertices of the feasible region, the equation of the hyperplane passing 

through them is 

4(Pi + P2 + P3 + Pa) = 1 (3.32) 
This hyperplane is not a part of the boundary of the feasible region because maximizing the left 
hand side of the above equation with respect to the product states shows that the point (17-, tt-, tt-, z) 



corresponding to the product state 



1 
1 

V 1 J 



1 
1 

V 1 / 



does not lie on the hyperplane. So the mentioned point is another vertex of the feasible region. 
Next, if we take the vertices {\, 0, 0, 0), (0, \, 0, 0), (0, 0, J, 0) and (-L, -L, i, I), the equation of the 
hyperplane passing through them is 

4(pi+P2+P3)+P4 = l- (3-33) 
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This is an exact boundary hyperplane of the feasible region. Therefore, any separable state must 
satisfy 4(pi + p 2 + ps) + Pa < 1. The existence of such hyperplane as a part of the boundary of the 
feasible region allows us to use the approach of the previous section to construct EWs. However, due 
to the lose of intuition in the four dimensional case, we can not deal with this problem analytically 
as well as in the 3 (g> 3 case and we have to invoke to the numerical evaluation. Let us consider the 
three vertices (0, |, 0, 0), (0, 0, |, 0) and ^, |) lying on the hyperplane (13.331) and the point 
(a,0,0,0) with a > |. The hyperplane passing through these points is 

v 1 1 

- + 4(p 2 + p 3 + (2 - — )p 4 = 1. (3.34) 

a 4a 

Numerical evaluation shows that when | < a < |, the hyperplane remains in contact, without 
intersecting, with the feasible region at the restricted plane passing through the vertices (0, 4, 0, 0), 
(0, 0, ~, 0) and (i, i, i, |). Thus, the inequality ^ + 4(p 2 +p 3 ) + (2 - ^)p 4 < 1 must be satisfied 
for any separable state. For a > |, the mentioned hyperplane does not remain in contact with the 
feasible region because it intersects the feasible region. In comparison with the 3 <g) 3 case, this 
situation, i.e. the variation of a in the mentioned range, is similar to the rotation of the plane (12.141) 
around the line segment which was discussed in the previous section. Permutation of the particles 
transforms the hyperplane (|3.34l) into the following one 

- + 4(p a + P2 ) + (2 - ±-)pt = 1. (3.35) 
a 4a 

In this case, numerical evaluation shows that, as for the hyperplane (13.341) . when | < a < | the 
hyperplane remains in contact with the feasible region at the restricted plane passing through the 
vertices (|, 0, 0, 0), (0, |, 0, 0) and (jq,jq,jqi |)- Hence, for any separable state we have ^ + 4(pi + 
p 2 ) + (2 — 4^)^4 < 1. F° r Qt > |, the hyperplane intersects the feasible region. If we take the vertices 
(|, 0, 0, 0), (0, 0, |, 0) and (y^, ^, |) and the point (0, a, 0, 0) with a > |, the passing hyperplane 
through them is 

^ + 4(p! + p 3 ) + (2 - -l)p 4 = 1 (3.36) 
a 4a 
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This hyperplane is invariant under the permutation of particles. Numerical results show that for all 
a > j, the hyperplane in (13.361) remains tangent to the feasible region, hence any separable state 
satisfies the inequality — + 4(p 1 + p 3 ) + (2 — -^)Pa < 1- The contact area is the restricted plane 
passing through the mentioned three vertices. Finally, the hyperplane passing through the other 
vertices (J, 0, 0, 0), (0, \, 0, 0), (0, 0, \, 0) and the point (0, 0, 0, a) is 

4(pi+p 2 +p 3 ) + - = l (3-37) 
a 

wherea G (— oo, 0) (J[l, oo). Fora = 0, we have the hyper-plane p 4 = as a part of the boundary of 
the feasible region because every separable state satisfies the p 4 > 0. When a varies in this range, the 
corresponding hyper-planes remain in contact with the feasible region at the restricted plane passing 
through the vertices (~, 0, 0, 0), (0, 4, 0, 0) and (0, 0, |, 0). Now we are ready to construct our EWs. 
The EW corresponding to the hyperplane (13.341) is 

W Q = h ® h - -O x - 4(0 2 + 3 ) - (2 - -1)04, (3.38) 
a Aa 

Permutating of the particles gives the EW 

W' a = h®h- -0 3 - 4(0! + 2 ) - (2 - -1)0 4 , (3.39) 
a 4a 

with | < a < |, which corresponds to the hyperplane (13.351) . We claim that the W a and W' a are 
nd-EWs. To prove this, we introduce the following state 

p = a x O x + a 2 2 + a 3 3 + a 4 4 (3.40) 

in which > 0,i = 1,2, 3, 4 and a\ + a 2 + a 3 + a 4 = 1. The PPT conditions for this state are 
a\a 3 > a\ and a 2 > a 4 . The expectation values of W a and W a with respect to the p are respectively 

Tr(W a p) = (oi - a 4 )(l - ^) (3.41) 

4a 



and 



rr«p) = (o3-a 4 )(l-^) (3.42) 
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As a special case of the state (13.401) . we consider the following state 

Q = J^—O, + j^—0 2 + ^—^0 3 + -r^—O, (3.43) 

13 + 7 13 + 7 13 + 7 13 + 7 

where < < 10 and 7 is non-negative. It is PPT for 1 < < 9 and 7 > 3 and not PPT for < 1, 
> 9 or 7 < 3. Taking the expectation values of W a and W a with q shows that q is PPT entangled 
state for 1 < < 3, 7 < < 9 and 7 > 3. For < < 1 and 9 < < 10, it is free entangled and 
for 3 < < 7 and 7 > 3 we can not say anything about its separability. The EW corresponding to 
the hyperplane (13.361 ) is 

< = h <8> h ~ -0 2 - 4(0! + 3 ) - (2 - -L)0 4 , (3.44) 
a 4a 

whose expectation value with the g shows that it is free entangled for 7 < 3 and does not detect it in 

the PPT regime. Finally, the EW corresponding to the hyper-plane (13.371) is 

WC = I 4 ® h ~ 4(Oi + 2 + 3 ) - -0 4 (3.45) 

a 

For a G (— 00, 0) |J[l,oo) it is always a positive operator and for a = 0, we have Wq = 4 
which is also a positive operator. To give a physical meaning for the EWs W a and W a / in detecting 
the entanglement of a state experimentally, from a measurement point of view, as in the previous 
section, these EWs should be decomposed into a sum of locally measurable operators. There exist 
following decompositions for W a and W a > as 

W a = ^ — -h ®h + 8 " X (A 2 ® A 2 - Ai <g) Ai + A 5 ® A 5 - A 4 ® A 4 + A 7 ® A 7 - A 6 ® A 6 
64a 32a 

3 

+A10 <g> A10 - A 9 <g> A 9 + A12 ® A12 - An <S> An + A14 <8> A u - A i3 ® A i3 ) + — — A 3 ® A 3 

32a 

16a + 5 x . 8a + 7, . y/Z{4a-\) f . . . 

+—^- R A 8 g> A 8 + — — A15 <8> Ai 5 + — ^— ^(A 3 g> A 8 - A 8 <g> A 3 ) 

96a 96a lb a 

\/2(4a-l) /x , , \/6(l-4a)^ , . 

+ ^— ^ - A 8 ® Ai 5 - Ai 5 ® A 8 + KA15 ® A 3 (3.46) 

12a 24a 

and 

W a > = 24 " - 3 / 4 ® £t + 8 " / (A 2 ® A 2 - Ai ® Ai + A 5 ® A 5 - A 4 ® A 4 + A 7 ® A 7 - A 6 ® A 6 
64a' 32a' 
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3 

+Ai <g> A10 - A 9 <g> A 9 + A12 ® A12 - An ® An + Au <8> A i4 - A i3 <g> A13) + ^—7^3 ® A 3 

62a 

16a' + 5, , 8a' + 7, , \/3(4a' - 1) , , 

+ nr ; A 8 <g> A 8 + — A15 ® A15 + ^— — - 3A 8 g> A 3 - A 3 <8) A 8 

96a' 96a' 48a' 

+ 1} (2Ai 5 ® A 8 - A 8 ® A 15 ) + V ^ ( 2 1 4 ~ (A3 ® Ai 5 ), (3.47) 

where A«s are the basis for the su(A) Lie algebra [20]. Therefore to measure the EWs W a and W a >, 
Alice and Bob need twenty local measurement settings. At the end of the paper, it should be noted 
that in order to characterize entanglement properties of other PPT states in 3®3 and 4®4 systems, and 
for higher dimensional bipartite one, several efforts have been made elsewhere(see e.g. ETll22ll23l ). 



4 Conclusions 

We have constructed new sets of EWs for the 3 <8> 3 systems analytically and for the 4 <g> 4 ones 
numerically. These witnesses are able to detect entanglement of some states in PPT regime. It 
has been shown that the existence of such EWs directly depends on the existence of a plane or a 
hyperplane as an exact boundary of the feasible region for the 3<g)3 and 4®4 systems respectively. The 
intersections of the plane or hyperplane with other boundaries of the feasible region are line segments 
or restricted planes respectively. As described, the rotation of the plane or hyperplane around the line 
segments or restricted planes in the allowed range leaves them in contact with, without intersecting, 
the related feasible region at the mentioned line segments or restricted planes. This situation, in fact, 
is due to the non-differentiability of the boundaries of the feasible region at the intersections. In this 
way, we obtain an infinite number of EWs. The approach can be generalized to the n ® n systems. To 
this aim, one can choose the operators O x — ^ J^tTo (^»* ® S)\ii)(ii\(I n <8> S' 1 "), 2 = \ Y^=o(I n ® 
S 2 )\ii)(ii\{I n ® S 2 t),..., O n _i = ^IXoVn ® S<> n -V)\ti){ii\(I n ® S^-W) and O n = as a 

basis for constructing EWs where S is the shift operator. Therefore, it is expected that there exists a 

hyperplane as an exact boundary of the feasible region of the form n(p 1 +p 2 +P3-\ \-p n ~i)+Pn = 1- 

The existence of such hyperplane motivates us to construct similar witnesses as constructed for the 
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3® 3 and 4® 4. If such witnesses are constructed, they will detect the entanglement of the generalized 
Horodecki state, i.e. p = YH=i a tOi, in the PPT regime. 



Pi 


P2 


^3 


Product state 


1 

3 


U 


n 

u 


v.rvJ 


u 


l 
3 


u 


{.)*{>) 


U 


U 


i 

3 


v.rvJ 


1 

9 


1 
9 


1 

3 




1 

48 


3 

16 


1 

4 




1 

192 


49 
192 


7 

48 




3 

64 


25 
192 


5 

16 


\\ ^ ) ® i ( ^ ) 

V vTo / \ VB / 


1 

12 


1 
8 


1 

3 


no product 


3 
16 


1 

48 


1 

4 


/ (J \ / Q \ 


49 
192 


1 
192 


7 

48 


3 f ^ ) ® i ( ^ J 


25 
192 


3 

64 


5 

16 


\ V6 / V ^ J 


1 

8 


1 

12 


1 

3 


no product 



Table 1 : The product states and coordinates of vertices for 3 3 EWs. 



References 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge 
University Press, Cambridge, 2000). 



Entanglement Witnesses 18 

[2] The Physics of Quantum Information: Quantum Cryptography, Quantum Teleportation and 
Quantum Computation, edited by D. Bouwmeester, A. Ekert, and A. Zeilinger (Springer, New 
York, 2000). 

[3] J. Preskill, The Theory of Quantum Information and Quantum Computation 

(California Inatitute of Technology, Pasadena, CA, 2000), 

http://www.theory.caltech.edu/poeole/preskill/ph229/. 

[4] A. Peres, Phys. Rew. Lett. 77, 1413 (1996). 

[5] P. Horodecki, Phys. Lett. A 232, 333 (1997). 

[6] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Rew. Lett. 80, 5239, (1998). 

[7] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 223, 1 (1996). 

[8] B.M. Terhal, Phys. Lett. A 271, 319 (2000); Linear Algebr. Appl. 323, 61 (2000). 

[9] S. L. Woronowicz, Rep. on Math. Phys. 10, 165 (1976). 
[10] S. Boyd, L. Vandenberghe, Convex Optimization (Cambridge University Press, 2004). 
[11] E.K.P. Chong, S.H. Z ak, An Introduction to Optimization (John Wiley, NY, 2001) 
[12] M. A. Jafarizadeh, M. Rezaee, S. K. A. Seyed Yagoobi, Phys. Rev. A 72, 062106 (2005). 
[13] M. A. Jafarizadeh, R. Sufiani, Phys. Rev. A 77, 012105 (2008). 

[14] M. A. Jafarizadeh, G. Najarbashi, Y. Akbari, H. Habibian, Eur. Phys. J. D 47, 233-255 (2008). 
[15] M. A. Jafarizadeh, Y. Akbari, N. Behzadi, Eur. Phys. J. D 47, 283-293 (2008). 
[16] M. Horodecki, P. Horodecki, R. Horodecki, Phys. Rev. Lett. 82, 1056 (1999) 
[17] A.C. Doherty, PA. Parrilo, F.M. Spedalieri, Phys. Rev. Lett. 88, 187904 (2002). 



Entanglement Witnesses 19 

[18] O. Giihne, P. Hyllus, International Journal of Theoretical Physics, Vol. 42, No. 5, May 2003. 

[19] G. Toth, O. Giihne, Phys. Rev. Lett. 94, 060501 (2005). 

[20] W. Pfeifer, The Lie Algebras su(N) An Introduction (Birkhauser Verlag, 2003). 

[21] H. P. Breuer, Phys. Rev. A. 71, 062330 (2005). 

[22] H. P. Breuer, J. Phys. A: Math. Gen. 38, 90199037 (2005). 

[23] H. P. Breuer, Phys. Rew. Lett. 97, 080501 (2006). 



